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Abstract 

We compute the influence of Coulomb effects on the charged kaon distribution at low 
momentum in high energy nucleus-nucleus collisions. This is accomplished by solving the 
Vlasov equation for the kaons in the presence of an expanding, highly charged, fireball. 
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The main goal of colliding heavy nuclei at high energy is to produce high density, 
high temperature matter, either hadronic or quark-gluon plasma. One way of studying 
the properties of high energy density matter is to compare the single-particle momen¬ 
tum distributions to those in a nucleon-nucleon collision at the same beam energy. The 
many-body medium should be most influential on the low momentum, or long deBroglie 
wavelength, part of the spectra and least influential on the high momentum, or short 
deBroglie wavelength, part. For example, there was a preliminary report of an enhance¬ 
ment of the charged kaon distributions in the low momentum region as compared to 
nucleon-nucleon collisions in an experiment at Brookhaven National Laboratory’s AGS 
[2]. Although this enhancement is now known to have been caused by an albedo source at 
the spectrometer’s collimator edge close to the beam trajectory 0, it begs the question 
of how the enormous electrical charge in a central gold-gold collision (Z = 158) modifies 
the low momentum spectra of charged hadrons. 

The most widely used approach to incorporate Coulomb effects into the description 
of the charged particle distributions is by means of the Gamow factor 0. 


G(v) 



( 1 ) 


Here rj = ZiZ 2 a/v ) v is the relative speed of the particle (charge Z{) with the source 
(charge Z 2 ), and ^(r = 0) is the relative Coulomb wave function evaluated at zero sepa¬ 
ration between the two (point) charges. In this approach, the charged particle distribution 
is given in terms of the neutral particle distribution and G by 
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where k(p) is the shift in the momentum due to the Coulomb kick. 

The original motivation for the introduction of eq. (i is the similarity to barrier 
penetration in fission and because this formula arises in a perturbative calculation to first 
order in Za for low momenta. But for heavy ion reactions, where the total charge Z is such 
that Za ~ 1, a first order perturbative approach to Coulomb effects is likely to fail and so 
will an approach that mimics such. Figure 1 shows the neutral, positive, and negative kaon 
momentum distributions in the fireball rest frame obtained from eq. (0) using a Maxwcll- 
Boltzmann distribution for neutrals and a non-expanding fireball with charge Z = 158 
and kaon temperature T = 120 MeV. The neutral distribution is normalized to unity. 
Generally the Gamow-modifed charge distributions do not preserve the normalization so 
the areas under those curves are not 1. Figure 2 shows the exact classical result obtained 
from the formula 
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where Vc(r) is the Coulomb potential felt by the kaon. The initial phase space distribution 
/ is taken to be a Maxwell-Boltzmann equilibrium distribution 

f r i v , to) oc exp{-[mu 2 /2 + V c (r)]/T} <d(R 0 - r) (4) 


which constrains the kaons to lie within a uniform density charged fireball of radius Rq = 8 
fm. The differences between the two figures are quite apparent and indicates the extent 
to which the Gamow-modihed neutral spectrum can be relied on. 

A more important physical effect which the Gamow factor (1) does not take into 
account at all is that in a high energy heavy ion collision the matter blows apart, and in 
fact it is the low energy detected charged particle which is left behind. 

Therefore it is desirable to look for an alternative approach to the description of 
Coulomb effects on charged particle spectra in high energy heavy ion collisions. In this 
paper we explore the influence of a time-dependent electric held produced by an expand¬ 
ing charged fireball on the low momentum spectrum of a test charge particle. To be 
suitable as a test particle the hadron must not be produced in great abundance and its 
nonelectromagnetic interaction with other hadrons at low relative momentum must not 
be very strong, otherwise it will influence the expansion of the fireball and the problem 
must be solved self -consistently. Kaons are good test particles at AGS type energies. For 
simplicity and definiteness we assume a thermal momentum distribution at the strong 
interaction freeze out and dynamically evolve this distribution by means of the Vlasov 
equation. The interaction of a single test particle with the rest of the matter is purely 
electromagnetic. Since the particles that produce the electric held are moving away from 
the collision region, the external held is time dependent. We describe the motion of non— 
relativistic test particles, thus ignoring interactions with the magnetic component of the 
electromagnetic held. We model the expanding hreball as a uniformly charged sphere, as 
viewed from the center-of-mass of the colliding nuclei, with a radius growing linearly with 
time. This expansion model is motivated by the similarity solution to a hydrodynamically 
expanding hreball [[| . This model turns out to be solvable analytically albeit in terms of 
Bessel functions. A serious comparison with experiment may require a simulation based 
on an event generator like ARC 0 or RQMD [0]. Finally, we remark that a description 
in terms of classical equations of motion for the test charge should be adequate for kaons 
with kinetic energy greater than about 10 MeV. At this energy the kaon DeBroglic wave¬ 
length is 2 fm which is small enough compared to expected gradients of temperature and 
density in the hreball. 

A uniformly charged sphere which has a radius i? 0 at time t 0l whose total charge is 
Ze, and whose radius R is increasing with time at a constant speed v s produces an electric 
potential 


V{r,t) 


Ze/Anr , r > R = v s t 

Ze(3R 2 — r 2 )/8nR 3 , r < R — v s t 


( 5 ) 


The hreball parameters are not independent but are related by Ro = v s t 0 . We hrst 
concentrate on the interior (r < v s t) potential. If / ± (r, v, £) represents the ±e charge 
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particle phase space distribution then, when ignoring particle collisions after decoupling, 
its dynamics is governed by Vlasov’s equation. In the interior region 


d d t s r d 

dt dr At 3 <9v 


/ ± (r.v,i) = 0, 


( 6 ) 


where 


t* = 


Ze 2 


7T mv 


( 7 ) 
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m is the kaon’s mass, and t > t 0 . Notice that c, the speed of light, does not enter into 
the definition of the characteristic time t s since we work in the non-relativistic limit. 

The solution to eq. (|j) is found by the method of characteristics. This involves solving 
the classical equations of motion and using the solutions to evolve the initial distribution 
/ ± (r, v, t 0 ) forward in time. 

In the interior the equations of motion are 


= v. 


or equivalently 


dr 

dt 

dv t s 

dt ~ ± 4f3 r ’ 
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-= ±—r . 
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For the initial distribution we take a Maxwell-Boltzmann equilibrium distribution at 
temperature T. 


/ ± (r, v, t 0 ) oc exp{ —[mu 2 /2 ± eV(r, t 0 )]/T} O(R 0 - r) 

Note that this constrains the kaons to lie within the initial sphere. 

The solution to the equations of motion for positive kaons is 


( 10 ) 


r(f) = A Vth (Jtjt) + BVtKi (Jti/t) 


and for negative kaons is 


r(t) = AVtJi + Bx/tiVi i^JtJPj 


( 11 ) 
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where the constant vectors A and B are determined by the initial conditions. The solution 
to eq. (|) incorporating the initial condition is 


/ ± ( r , v, t) = N exp < —— 
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The normalization constant J\f is chosen so that the distribution is normalized to unity 
at the initial time to- 

J d 3 r d 3 v f ± (r,v,t 0 ) = 1 (14) 

Of course, this ensures that the distribution is normalized to unity for all later times, even 
when account is taken of the fact that the solution differs outside the expanding sphere. 
For negative kaons the dimensionless functions of time Cj(t) are given by 
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For positive kaons the corresponding Cj(i) are obtained by the substitution t s 
the right-hand side of eqs. ( |15|) and the identities 


(15) 
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Generally the asymptotic momentum distribution is obtained by integrating over po¬ 
sition and taking the limit t —> 00. Once a kaon has crossed the surface of the expanding 
sphere it will see a time-independent potential, hence its energy will thereafter be con¬ 
served. Since a positive kaon will see a positive potential, its kinetic energy will not 
decrease, and neither will its velocity. Once a positive kaon has crossed the surface it will 
not come back. This will not be true for negative kaons; they will be discussed later. Due 
to spherical symmetry, the distribution is a function only of the magnitude of the three 
dimensional momentum vector. The coordinate integration can be performed analytically 
with the result for positive kaons being 
d 3 N 
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where $(x) is the error function and 
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where the Bessel functions are all evaluated at yt s /to- Taking Z = 158, Ro = 8 fm, v s /c 
= 0.4, and T = 120 MeV, results in t 0 = 20 fm/c, t s = 28.7 fm/c, pc = 1390 MeV/c, x 
= 1.19, and S = 0.330. It is important to remember that eq. (0) is valid only so long as 
the test particle (positive kaon) remains within the expanding sphere. This will be the 
case if its asymptotic speed is less than v s . 

To compute the distribution of positive kaons for asymptotic speeds greater than v s 
is straightforward due to the fact that once it has crossed the boundary of the expanding 
sphere, it sees a time independent potential so that its energy is thereafter conserved. 
The flux of particles crossing the surface at time t is 

J d 3 rv ■ r/ + (r, v,t)\ r=Vst 


which is to be evaluated at r = v s t. The energy distribution is obtained by multiplying 
by the surface area, inserting an energy conserving 6 function, and integrating over time. 

POO p p 

dN + /dE = / dtv 2 t 2 / dfl r / d 3 v v ■ r/ + (r, v,t)\ r=Vst 

Jt 0 J J 

x 5(E — p 1 /2m — Ze 2 /4:7Tv s t) (19) 

The asymptotic momentum distribution is readily obtained from the energy distribution. 

Figure 3 shows the K + distributions for expanding fireballs with initial radius 8 fm, 
kaon temperature 120 MeV, and expansion velocities v s /c = 0.2, 0.4 and 0.6. For compar¬ 
ison the dotted curve shows the asymptotic distribution obtained from a static fireball. 
There is a discontinuity in the slope of the distributions at E — \mv 2 . We have checked 
that the areas under all curves are 1. The significant physics here is that the distribution 
remains finite at zero energy. The reason is that the charged fireball blows apart and 
leaves a finite number of low energy kaons behind. This cannot happen with a static 
fireball where a classical turning point exists. In the limit that the fireball expansion 
velocity is large compared to the thermal velocity of the kaons, the kaon distribution is 
hardly changed from the pure exponential form. The fireball expands so fast that it leaves 
the kaon distribution frozen to its initial shape. In the limit that the expansion velocity 
goes to zero, the static fireball result is recovered. 

Now we turn to negative kaons. A negatively charged kaon can cross the fireball’s 
surface if it has a sufficiently large radial velocity. After crossing the surface, it will 
slow down due to the attractive potential. It may even slow enough so that the surface 
passes it by! So the possibilities are three-fold: it always remains within the expanding 
sphere, it crosses the surface once and never comes back, or it crosses the surface twice so 
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that asymptotically it remains within the expanding sphere. Kaons with an asymptotic 
energy greater than | mv 2 can be dealt with by the flux formula (|19D. Kaons with an 
asymptotic energy less than this cannot be handled analytically. This is ultimately due 
to the fact that the Coulomb trajectories cannot be expressed in terms of elementary 
functions of the time and so it is not possible to construct the phase space distribution 
outside the expanding sphere in a useful form. To handle this situation we have written 
some computer programs which calculate the trajectory of a kaon given its initial phase 
space position. One program was written in turbopascal with a screen display in real time 
which graphically shows the behaviors sketched above. A second program was written 
in fortran which connects the initial phase space point to the final energy (usually the 
asymptotic energy was determined at lOOfo)- 

Figure 4 shows the K~ distributions corresponding to the same fireball configurations 
as in figure 3. The solid curves show the result of the flux calculation for asymptotic 
energies greater than ^mv 2 . The histograms are the result of numerically computing 
the kaon trajectories for 720,000 initial phase space points. The initial radial position 
was incremented in 0.1 fm bins, the initial speed was incremented in 0.01c bins, and the 
angle between the initial position and velocity vectors was incremented in 2° bins. (Hence 
80 x 100 x 90 = 720, 000.) The distribution is enhanced at small energies relative to the 
neutral kaons and to the positive kaons, but not as much as for a static fireball, shown 
in the figure by dotted curves. Another distinction is that there are no classically bound 
J7 _ , unlike the case of a static fireball. 

The kaon data which originally attracted attention was for rather forward rapidities 
of 2.1 to 2.5 in contrast to the rapidity of the nucleus-nucleus center-of-mass frame of 
about 1.5. Plotted in the nucleus-nucleus center-of-mass frame the kinetic energies go no 
lower than 80 MeV for K~ and 70 MeV for K + . It is problematic to compare against 
theoretically calculated spectra for kinetic energies greater than about 100 MeV when a 
nonrelativistic approximation has been made. In addition, the selection criterion was for 
the 10% most central collisions, corresponding to the 30% most central impact parameters. 
These impact parameters are really too large to focus on truly central collisions of equal 
mass nuclei. 

In conclusion we have shown the importance of the Coulomb force on charged kaons 
in central collisions of large nuclei at high energies. We have shown how the expansion 
of the fireball affects the kaon distributions in a quantitative way. Comparison to data 
at kinetic energies greater than about 100 MeV probably requires a relativistic treatment 
which we are now working on. By varying the initial phase space distribution of the kaons 
by invoking some component of flow, for example, and comparing with data should yield 
quantitative information on the dynamics of kaons and heavy ion collisions. 
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Figure Captions 

Fig. 1: Gamow factor momentum distributions d 3 N/d 3 p as a function of kinetic energy 
E for a non-expanding fireball with T = 120 MeV and Z = 158. Also shown is the 
corresponding neutral distribution. 

Fig. 2: Momentum distributions d 3 N/d 3 p as a function of kinetic energy E for a non¬ 
expanding fireball with Rq = 8 fm, T = 120 MeV, Z = 158. The sudden drop in the K + 
distribution signals the proximity of E to E min = Ze 2 /47i\R 0 , the classical turning point. 
Compare this to figure 1. 

Fig. 3: K + momentum distributions d 3 N/d 3 p as a function of kinetic energy E for a 
uniformly expanding fireball with R 0 = 8 fm, T = 120 MeV, Z = 158 and various values 
of the expansion velocity v s indicated with solid lines. Also shown is the corresponding 
result for a static fireball v s = 0 indicated with dotted lines. 

Fig. 4: K~ momentum distributions d 3 N/d 3 p as a function of kinetic energy E for a 
uniformly expanding fireball with Rq = 8 fm, T = 120 MeV, Z = 158 and various values 
of the expansion velocity v s . The solid lines starting at E s = are obtained from 

the flux across the fireball’s surface. The histograms are obtained by solving the classical 
equations of motion numerically for 720,000 starting points in phase space. Also shown 
is the corresponding result for a static fireball v s = 0 indicated with dotted lines. 
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